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ABSTRACT 

In this paper we answer a simple question: can a misaligned circumbinary planet induce 
Kozai-Lidov cycles on an inner stellar binary? We use known analytic equations to analyse 
the behaviour of the Kozai-Lidov effect as the outer mass is made small. We demonstrate a 
significant departure from the traditional symmetry, critical angles and amplitude of the effect. 
Aside from massive planets on near-polar orbits, circumbinary planetary systems are devoid 
of Kozai-Lidov cycles. This has positive implications for the existence of highly misaligned 
circumbinary planets: an observationally unexplored and theoretically important parameter 
space. 

Key words: binaries: close - astrometry and celestial mechanics: celestial mechanics - plan¬ 
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1 INTRODUCTION 

In 1962 there were two independent papers published on an effect 
seen in the gravitational three-body problem. Lidov (1962) inves¬ 
tigated the effect the Moon has on artificial satellites orbiting the 
Earth. Kozai (1962) looked at asteroids orbiting the Sun under the 
influence of Jupiter. Both systems were qualitatively the same: an 
inner restricted three-body problem with initially circular orbits and 
all of the orbital angular momentum confined to the outer orbit. 
(Fig. la). It was discovered that when the two orbital planes were 
significantly misaligned, between 39° and 141°, there was a high- 
amplitude modulation in the eccentricity of the inner orbit (in their 
examples the satellite and the asteroid) and the mutual inclination. 
Contrastingly, the outer eccentricity was seen to remain constant. 
Overall, this is known as a Kozai-Lidov (K-L) cycle or effect. 

The K-L effect was later generalised to the case of three bod¬ 
ies of comparable mass (Harrington 1968; Lidov & Ziglin 1976) for 
the application to triple star systems (Fig. lb). It has subsequently 
been implicated, for example, in the production of tight stellar bi¬ 
naries (period < 7 d). In this scenario the K-L cycle induced on 
the inner binary leads to close encounters between the two stars, at 
which point tidal friction dissipates orbital energy and shrinks the 
inner orbit (e.g. Mazeh & Shaham 1979; Eggleton & Kisseleva- 
Eggleton 2006; Fabrycky & Tremaine 2007). 

The theory has been extended to eccentric outer orbits, in 
which case higher-order effects cause the system to be chaotic, pos¬ 
sibly inducing flips in the inner orbit orientation. This has been ap¬ 
plied to multi-planet systems (Fig. Ic) by Naoz et al. (201 1) in order 
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to explain why a surprisingly large fraction of hot-Jupiters are mis¬ 
aligned or even retrograde (Hebrard et al. 2008; Schlaufman 2010; 
Triaud et al. 2010). 

The similarity between the three scenarios in Fig. la - c is that 
the system’s orbital angular momentum is largely confined to the 
outer orbit. The timescale, amplitude and limits of K-L cycles do 
vary between the different configurations, but for initially circular 
orbits there is qualitatively not a radical departure from the original 
work of Lidov (1962) and Kozai (1962). 

In this paper we explore a different type of three-body sys¬ 
tem: circumbinary planets (Fig. Id). The discovery of eight tran¬ 
siting systems by the Kepler mission (Doyle et al. 2011; Welsh et 
al. 2014) has given this class of planet widespread scientific expo¬ 
sure and validity. The discoveries so far have been limited to copla- 
nar systems (within ~ 4°), so one might question the relevance of 
K-L cycles in this scenario. However, Martin & Triaud (2014) ar¬ 
gued that this narrow distribution is the result of a strong detection 
bias imposed by the requirement of a consecutive transit sequence^. 
There are also theoretical arguments for the existence of misaligned 
circumbinary planets, as a result of planet-planet scattering (e.g. 
Chatterjee et al. 2008 in the context of single stars), a misaligned 
disc (e.g. 99 Herculis, Kennedy et al. 2012, KH 15D, Winn et al. 
2004) or an interaction with an outer tertiary star (Munoz & Lai 
2015; Martin et al. 2015; Hamers et al. 2015). 

In a circumbinary planetary system the outer orbital angular 
momentum becomes vanishingly insignificant with a decreasing 


^ Circumbinary planets that are misaligned by > 10° frequently miss tran¬ 
sits, creating a sparse transit sequence which is harder to identify. 
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Figure 1. Zoo of different three-body configurations, all with orbits misaligned by A/. The orbits are drawn with respect to the centre of mass of the inner 
orbit. Case (d) is qualitatively different to cases (a) - (c) because in (d) the outer orbit no-longer dominates the total angular momentum of the system. 


outer mass. This is qualitatively different to the other systems in 
Fig. 1 and it means that the classic studies of Lidov (1962) and 
Kozai (1962) are no-longer entirely applicable. Furthermore, since 
the K-L effect only modulates the eccentricity of the inner orbit, 
not the outer, it is not possible for an inner binary to induce a K- 
L cycle on a circumbinary planet. Naoz et al. (2013b) and Liu et 
al. (2015) derived equations for the amplitude and limits of the K- 
L effect for three arbitrary mass bodies, but did not apply them to 
circumbinary planets. We are motivated to demonstrate the clear 
and significant dijference between circumbinary planets and other 
three-body systems, and to clear up any possible misconceptions. 

The plan of the paper is as follows: In Sect. 2 we present equa¬ 
tions for K-L theory, analyse the different limiting cases and con¬ 
sider the competing secular effect of general relativistic precession, 
which may act to suppress K-L cycles. We then apply the theory in 
Sect. 3 to determine what parameter space of circumbinary planets 
may be able to induce K-L cycles before concluding in Sect. 4. 


2 THEORY 

2.1 Equations for the induced eccentricity and critical 
mutual inclination 

Consider a hierarchical triple system like those shown in Fig. 1, 
which we model as two Keplerian binary orbits. The three bod¬ 
ies have masses Mi (navy blue in Fig. 1), M 2 (blue) and M 3 (light 
blue). Both orbits are defined by osculating orbital elements for the 
period, P, semi-major axis, a, eccentricity, e, argument of periapse, 
cx), inclination I and longitude of the ascending node, Q. We use the 
subscripts “in” and “out” to denote each orbit. The mutual inclina¬ 
tion between the two orbits. A/, is calculated by 


cos A/ = sin /in sin /out cos AQ -h cos /in cos 4^, (1) 

where AQ = - flout- The orbital angular momenta of the inner 

binary and outer binaries are 
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where 


MiMo I - 

Ml + M 2 


and 


(Ml -|- M2)M3 / -—-—-— 

Tout = — - — - — ^^G{Mi + M 2 + M2)aoy 

Ml + M2 + M2, 

and G is the gravitational constant. 


( 3 ) 


( 4 ) 


We analyse the three-body orbital evolution according to the 
quadrupole approximation of the Hamiltonian. In this case is no 
change to ain, ^out or ^out Under certain conditions there may be a 
significant variation in the ^in and A/: a K-L cycle. Outside of K-L 
cycles ^in is constant. Additionally, both orbits will experience an 
apsidal and nodal precision (variations in tu and fl, respectively). 

The presence or absence of K-L cycles has significant impli¬ 
cations for the stability of a three body system. For circular and 
coplanar triple systems there exists a rule of thumb that ^out ^ 3ain 
for stable orbits (Dvorak 1986; Holman & Wiegert 1999; Mardling 
& Aarseth 2001). Around circular binaries, the stability limit gen¬ 
erally moves inwards as Al increases and is only a weak function 
of the stellar mass ratio (Doolin & Blundell 2011). If the inner bi¬ 
nary is eccentric, for example during a K-L cycle, then this stability 
limit is pushed outwards. An eccentric binary also makes the stabil¬ 
ity limit a more complex function of Al and the mass ratio (Doolin 
& Blundell 2011; Li et al. 2014). Seven of the known eight cir¬ 
cumbinary systems have ^in < 0.2, with Kepler-34 being the only 
highly eccentric case (^in = 0.5, Welsh et al. 2012). 

Naoz et al. (2013b) and Liu et al. (2015) derived a relation 
for the maximum inner eccentricity obtained, ^in,max, as a function 
of the starting mutual inclination, A/q, the ratio of orbital angular 
momenta and the outer eccentricity. 
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where it is assumed that the inner binary is initially circular. So¬ 
lutions for ^in,max Only exist in a “K-L active” region, which is 
bounded by lower and upper limits on the initial A/, which we call 
A/lower and A/upper - 

The limiting mutual inclination for K-L cycles to occur is cal¬ 
culated by setting ^in,max = 0 in Eq. 5, leaving a quadratic 


5 cos All 
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^ In the octupole level Hamiltonian there may be some variation in ^out 
but generally this is small. An interesting exception however was found by 
Li et al. (2014), who showed that an eccentric inner binary (^in > 0.4) may 
induce moderate eccentricity variations (up to ^out ~ 0-3) on a massless 
outer body if it has a near-polar orbit. 
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Figure 2. Maximum inner eccentricity obtained during K-L cycles according to Eq. 5, as a function of fl!out/«in and A/q, for six different values of the outer 
mass M 3 . The inner masses are set at M\ = IM© and M 2 = O.SMq. The initial eccentricities are both zero (and ^out is constant). The black regions correspond 
to where K-L cycles do not occur. These are separated from the K-L active regimes by a lower limit of A/q (white dashed line on the left-hand limit from Eq. 7) 
and an upper limit of A/q (green and/or pink dashed lines on the right-hand limit from Eq. 8 and/or Eq. 9, respectively). The light blue horizontal lines are 
upper limits on ao^xtlcim for K-L cycles to occur without being quenched by GR precession, for = 0.0655 AU (Pin = 5 d) as the lower line and a\n = 0.3041 
AU (Pin = 50 d) as the upper line. The light blue star symbol symbol in the top right subplot corresponds to the set of N-body simulations in Pig. 3. 


From solving the quadratic the lower limit of A/q is defined as 


cos A/iower ~ 


1 

To 


Lin 1 

a/I-'ou. 



The upper limit is simply the other solution to the quadratic, 


cos A/upper ~ 


1 

To 


Lin 1 


for 


An 1 


< 2 , 



( 8 ) 


but with the caveat that it is only defined when the outer orbit con¬ 
tains the majority of the angular momentum. Otherwise one must 
invoke a different upper limit, taken from Lidov & Ziglin (1976): 


cos A/upper = 




for 


Lin 1 


> 2 . 


(9) 


This second upper limit is not inherently respected by Eq. 5 when 
calculating ^in,max so it must be imposed manually. 

Liu et al. (2015) compared numerical simulations with the an¬ 
alytic calculations using Eq. 5, showing a very close match, but 
they only used Eq. 8 as an upper limit on A/q, not Eq. 9. However 


in their tests the outer orbit always had a significant portion of the 
angular momentum, and hence Eq. 9 was never applicable. 


2.2 Limiting cases 

Kozai-Lidov cycles were originally derived in the limit of the inner 
restricted three-body problem, in which case M 2 0 and hence 
Lin/Lout ^ 0 - In this case cos A/iower ^ V3/5, and for the upper 
limit we use Eq. 8 to get cos A/upper ^ - V3/5. This recovers the 
classic limiting mutual inclinations for K-L cycles: A/iower = 39.23° 
and A/upper = 140.79°. The calculation for ein,max reduces from Eq. 5 
to the simple well-known formula 


which is symmetric around A/q = 90°. 

The opposite limiting case is the outer restricted three-body 
problem, where M 3 ^ 0 , and hence Lin/Lout ^ 00 . In this limit 
cos A/lower ^ 0 and hence A/iower ^ 90°. This means that the pro¬ 
grade K-L active region disappears. Eor the upper limit we use Eq. 9 
to see that cos A/upper ^ 0 and hence A/upper ^ 90°, so the retro¬ 
grade region disappears too. We recover the trivial limit of the outer 
restricted three-body problem, where there can be no K-L cycles 
since the binary feels no influence from the orbiting body, matching 
the work of Pilat-Lohinger et al. (2003); Earago & Laskar (2010); 
Doolin & Blundell (2011). The only dynamical effect on the cir¬ 
cumbinary planet’s orbit is an apsidal and nodal precession, at a 
rate calculated by Earago & Laskar (2010). 
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Figure 3. Numerical simulations of K-L cycles, showing the evolution of 
A/ (top panel), ein (middle) and ^out (bottom). The parameters are Am = 0.5 
AU, <3out = 5 AU, Ml = IMq, M 2 = 0.5Mq, M 3 = 0.05M© and ^out and ein 
initially set at zero. Three different values of AIq were tested: 157° (light 
blue), 158° (blue) and 159° (navy blue). All other orbital angles were set 
to zero. The top panel includes a zoomed inset to show the small difference 
in starting conditions. The numerical values of ein,max are 0.9894 (A/q = 
157°), 0.9899 (A/o = 158°) and 1.63 x lO"'* (A/q = 159°), and hence in 
the last case there is no K-L cycle. The analytic predictions for ein,max are 
correct to within 0.04% error. 


2.3 Competing secular timescales 

Even if a circumbinary system is said to be in a K-L active regime 
according to Sect. 2, it is still possible that K-L cycles will not 
occur. This is because general relativity (GR) induces a competing 
precession on the inner binary, and out of the two effects the slower 
one will be suppressed^. The K-L timescale is 


= 1 /fAif " I M1 + M2 
sUin/ y GMf 


(11) 


which is taken from Fabrycky & Tremaine (2007) but has been 
converted to be a function of semi-major axes. The GR timescale is 
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where c is the speed of light (Fabrycky 2010). 

An approximate limit for the suppression of K-L cycles is 


^ Tidal and rotational bulges also induce a precession on the binary’s orbit 
but these effects are generally only significant for very tight binaries. 


Tk-l 

tgr 


( 13 ) 


This criterion assumes a sharp cut-off of K-L cycles, whereas in 
reality there is likely a smooth transition between a K-L dominated 
regime and a GR dominated regime. We are also ignoring the possi¬ 
bility of resonances between the two secular effects, which may in 
fact boost the modulation of ^in (Naoz et al. 2013a). Nevertheless, 
this criterion is deemed sufficient for the purposes of the simple 
analysis presented in this paper. Evaluating Eq. 13 using Eqs. 11 
and 12 yields 


^out 


M3 

ZZin 

lim 

(Ml + M2f 2G 


(14) 


This provides an upper limit on the tightness of circumbinary sys¬ 
tems for which K-L cycles may be excited, lest they be quenched 
by general relativistic precession. 


3 PARAMETER SPACE EXPLORATION 

We explore how the Kozai-Lidov effect behaves in an example 
triple system consisting of an inner stellar binary with Mi = IM© 
and M 2 = 0.5M© and an outer body with M 3 ranging from a star 
down to a planet. Different values of Mi and M 2 would change 
the limits and amplitude of the K-L effect, but in our example we 
use the mean stellar masses of the known circumbinary systems 
from Kepler, in order to make it as representative as possible. From 
Sect. 2.1 the amplitude and limits of K-L cycles are functions of the 
ratio of angular momentum, and consequently in our tests we only 
need to probe different values of the ratio ^out/^^m, not the individ¬ 
ual values. We start ei^ and ^out at 0, but ei^ may rise significantly 
during a K-L cycle. 

In Fig. 2 we calculate ein,max from Eq. 5 as a function of aout/^^m 
and A/o, for six different outer body masses: IM©, O.IM©, 0.05M©, 
O.OIM©, 0.005M© and O.OOIM©. The smallest value of aout/^^m was 
3, corresponding to the rough stability limit. The lower limit (left- 
hand limit) of the K-L active region is demarcated in a white dashed 
curve (Eq. 7). The two different upper limits (right-hand limits), 
chosen according to the value of M 3 and aout/^in, are depicted as 
green (Eq. 8 ) and pink (Eq. 9) dashed curves. 

The relative timescales of K-L cycles and GR precession set 
an upper limit in aout/^in on the K-L active region, as a function 
of ain (Eq. 14). In Fig. 2 two light blue horizontal lines denote the 
limiting ratio aout/^^in for ai^ - 0.0655 AU (Pin = 5 d, lower line) 
and tZin = 0.3041 AU (Pin = 50 d, upper line). These two example 
values of ai^ are not completely arbitrary as they cover the range of 
binaries known to host circumbinary planets (Welsh et al. 2014). 

For M 3 = IM© we have a triple star system like in Fig. la, 
which is well-modelled by a classic K-L regime with limits at 
39° and 141°. Here ^in,max is solely a function of A/q and is well- 
approximated by Eq. 10. The GR timescale only impinges upon 
the widest systems within this parameter space. As M 3 decreases to 
O.IM©, an M-dwarf, the K-L active regime is relatively unchanged 
for wide systems but shifts towards the retrograde region for tight 
systems. This asymmetry was noted by Liu et al. (2015). 

Within the brown-dwarf regime (M 3 = 0.05M©) there is 
a significant change in the K-L behaviour for tight triple sys¬ 
tems (tZout/^in < 15). The high-amplitude eccentricity modulations 
(^in,max ^ 0-9) are confined to retrograde orbits. Also the second 
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upper limit (Eq. 9) becomes applicable, causing the upper limit of 
the K-L active region to “turn around” and move back towards 90°. 

Within the planetary regime between roughly 1 and 10 Mjup, 
we see that K-L cycles become restricted to a narrow slither of 
angles either side of a polar A/q = 90° orbit. The amplitude also 
decreases as the outer body becomes smaller. Furthermore the K- 
L timescale is also so slow that GR precession suppresses it in all 
but the closest systems. For reference, the observed circumbinary 
planets have been found predominantly near critical stability limit 
(flout ~ 3flin), with low eccentricities (^out < 0 . 2 ), have mass M 3 < 
IMjup and are nearly coplanar (A/ < 4°) (Welsh et al. 2014). 

One interesting property is that the upper limit of K-L cycles 
calculated by Eq. 9 (pink dashed curve in Fig. 2) imposes a sharp 
transition between K-L active and inactive regions. As an example, 
for M 3 = 0.05Mq and flout/^^in = 10 there is a sharp transition at 
A/o = 159° (marked with a light blue star symbol in the top right 
subplot in Fig. 2). We demonstrate this with an N-body numerical 
solution of the exact equations of motion ^ in Fig. 3, showing the 
time-evolution of the osculating orbital elements A/, ^in and ^out 
for three different values of A/q: 157° (light blue), 158° (blue) and 
159° (navy blue). The N-body code does not include GR, which is 
reasonable since the K-L timescale is significantly faster in this ex¬ 
ample. There is a sharp transition between K-L active and inactive 
regions at A/q = 159°. The maximum numerical value matches 
the analytic prediction to within 0.04% error. 

Figure. 3 shows that ^out remains close to zero but not exactly. 
For binaries undergoing a K-L cycle the outer eccentricity experi¬ 
ences low-amplitude periodic rises, in phase with the K-L cycle of 
the inner binary. This is not a K-L cycle being induced on the outer 
body itself, but rather a response to the high-amplitude eccentricity 
modulation of the orbited binary. 

Relaxing the requirement of a circular outer orbit moves the 
sharp upper limit towards 90°. In the Fig. 3 example with ^out = 0.5 
there is a shift in A/upper from 159° to 144°, according to Eq. 9, 
which we verified in an N-body simulation. Like before, the simu¬ 
lation showed little variation in ^out over time. 


4 CONCLUSION 

The near-complete absence of potentially destabilising Kozai- 
Lidov cycles in the context of circumbinary planets has positive 
ramifications for the existence of misaligned systems. It frees us 
from the narrow confines of coplanarity and broadens the param¬ 
eter space for potential discoveries. Based on studies of the abun¬ 
dance of circumbinary planets, if a presently-hidden population of 
misaligned systems were to exist, this would imply that planets are 
surprisingly more abundant around two stars than one (Armstrong 
et al. 2014; Martin & Triaud 2014). 

There are also theoretical implications in the context of close 
binary formation by a combination of K-L cycles and tidal friction. 
Circumbinary planets cannot activate this mechanism. Circumbi¬ 
nary brown dwarfs may but only in close, retrograde orbits. 

Observational evidence of such misaligned systems is 
presently lacking, but fortunately there exists several methods for 
their discovery: eclipse timing variations (Borkovits et al. 2011), 
astrometry (Sahlmann et al. 2015) and sparse transits on both 
eclipsing and non-eclipsing binaries (Martin & Triaud 2014, 2015). 

^ REBOUND (Rein & Liu 2012, http://github.com/hannorein/rebound), 
with a 14-15th-order adaptive time step integrator (Rein & Spiegel 2015). 
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